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We study superconductivity in doped solid picene (C22H14) with linear response calculations of 
the phonon spectrum and electron-phonon (ep) interaction. We show that the coupling of the high- 
energy C bond-stretching phonons to the tt molecular orbitals for a doping of ~3 electrons per picene 
molecule is sufficiently strong to reproduce the experimental Tc of 18 K within Migdal Eliashberg 
theory. For hole doping, we predict a similar coupling leading to a maximum Tc of 6 K. However, we 
argue that, due to its molecular nature, picene may belong to the same class of strongly correlated 
ep superconductors as fuUerides. We propose several experimental tests for this hypothesis and 
suggest that intercalated hydrocarbons with different arrangements and numbers of benzene rings 
may be used to study the interplay between ep interaction and strong electronic correlations in the 
highly non-adiabatic limit. 

PACS numbers: 31.15.A-, 74.70. Kn, 63.20.kd 
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The field of superconductivity has witnessed several 
important discoveries in the last ten years that, together 
with new synthesis and manipulation techniques, have 
enormously advanced our understanding of this phe- 
nomenon. These include MgB2 with Tc of 39 K,[lj boron- 
doped diamond with Tc as high as 11.5 K,[2j intercalated 
graphite[l| with a Tc of 14.5 under pressure, Q and iron 
based superconductors Q with Tc of up to 56 K.0 How- 
ever, in order to be able to design superconductors with 
desired properties, additional efforts are needed to un- 
derstand and control the material parameters relevant 
for superconductivity in different classes of compounds. 

Carbon based compounds offer many possibilities in 
this regard because their structure can be easily ma- 
nipulated at the atomic level and they generally posses 
high phonon frequencies that have substantial electron- 
phonon (ep) coupling. Besides many conventional 
phonon-mediated superconductors, this class also com- 
prises alkali- intercalated fuUerides (A3 Ceo) with Tc of up 
to 40 K.0, Q These 7r-bonded molecular solids have a 
rich phase diagram determined by a non-trivial interplay 
between strong local ep interactions and on-site Coulomb 
correlations in a highly non-adiabatic regime. 0- Hi The 
physics of this strongly correlated ep superconductors 
is of great fundamental interest. Therefore, identifying 
other compounds in this class is highly desirable. For 
this reason, the discovery of superconductivity in K and 
Rb doped picene with Tc of up to 18 K by Mitsuhashi et 
at last year[l3| is of great significance. 

In this paper, we use linear response calculations of 
the phonon spectra and ep interaction 13| within the 



rigid band approximation (RBA) to show that alkali- 
doped picene could be another case of strongly corre- 
lated ep superconductors. Our results complement pre- 
vious theoretical and experimental studies of this com- 
pounds, which have focused on electronic properties and 
correlations. [iB-lHl We obtain a value of the ep coupling 



constant that is large enough to explain superconductiv- 
ity with Tc = 18 K in Migdal-Eliashberg theory. How- 
ever, since the bandwidths of the conduction and valence 
bands in solid picene are relatively small (^ 0.3-1 eV), 
we argue that a local approach that includes both the ep 
coupling and Coulomb correlations on an equal footing 
may be more appropriate and also provide parameters 
that are relevant to these type of models. Finally, we 
discuss possible experimental tests that could confirm or 
exclude the scenario we propose. 

Picene (C22H14) is a polycyclic aromatic hydrocar- 
bon formed by five benzene rings juxtaposed in an arm- 
chair arrangement. Solid picene is an insulator, with 
a gap of ~3.3 eV, which can form in either monoclinic 
or orthorhombic structure. [2^ Since the superconducting 
samples of Ref. |l2| have a monoclinic P2i structure, we 



employed experimental lattice parameters for undoped 
picene in this variant with a = 8.480,6 = 6.154, c ~ 
13.515 A, /3 = 90.46°. [H] We relaxed the internal coordi- 
nates so that every component of force was less than 10~^ 
Ry/au. This changed the C-C and C-H bond lengths by 
less than 2.5% and was needed in order to have a dynam- 
ically stable structure with real phonon frequencies. [2^ 
The unit cell contains two picene molecules placed in a 
herringbone arrangement in the xy plane (shown in the 
inset of Fig. [T]) which are then stacked along z direction. 
The shortest intermolecular C-C distance is 3.5 A within 
the xy plane and 3.9 A along z (i.e. along the length 
of the picene molecule). These intermolecular distances 
are more than twice the typical intramolecular nearest- 
neighbour C-C distance of 1.4 A. In the doped super- 
conducting crystals, the a axis expands while the b and 
c axes shrink. It is inferred from this that the dopants 
most likely occupy intra-layer positions since inter-layer 
intercalation would have caused the c axis to expand. 
However, the detailed experimental structure of doped 
picene is not available. For this reason, we simulated 
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FIG. 1: (Color online) Calculated LDA electronic structure 
of solid picene. In green (light gray) and blue (dark gray) 
we show the two HOMO and the four LUMO* derived bands, 
respectively, which are accessible by hole or electron doping, n 
is the maximum number of electrons per picene molecule that 
can be doped. The two insets show the unit cell of solid picene 
in the xy plane and its Brillouin Zone (BZ), respectively. 



22| is shown in the main panel of Fig. [TJ The 



zero of the energy in the figure is chosen as the mid- 
dle of the calculated '^2.4 eV gap that separates the the 
bonding valence bands from the anti-bonding conduction 
bands in undoped picene. All bands in the energy range 
shown in the figure have mostly C Pz character. They 
are lumped into small subsets with very small in- and 
out-of-plane dispersion reflecting the molecular nature of 
the solid. This is similar to the electronic structure of 
fuUerenes whose conduction bands also have C Pz char- 
acter and are very narrow. Q In superconducting picene 
(K/Rb3C22Hi4) the ~3 electrons donated from the alkali 
atoms to each molecule populate the 4 bands immedi- 
ately above the gap. The Fermi surface (not shown) is 
large without any significant two-dimensional nature and 
prominent nesting features. [l7i [21I These bands have a 
total bandwidth of ~ 0.3 eV [23| and derive from a 
combination of the two lowest unoccupied molecular or- 
bitals (LUMO and LUMO-l-1), which have a small energy 
separation of ^66 meV in the molecule. [l7j The energy 
separation between the two highest occupied molecular 
orbitals (HOMO and HOMO-1) is instead higher (- 170 
meV) so that the four corresponding bands in the solid 
are grouped in two subsets of two bands separated by a 
gap of ^ 0.1 eV. [l3| In the following, we will study within 
RB A the effect of doping electrons or holes into the states 
immediately above (LUMO*) and below (HOMO) the 
gap, which are colored in blue and green, respectively, in 
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the effect of doping using the rigid band approximation 
(RBA) and did not include cxplicitely the dopants in the 
calculations. The role of the dopant, which should be 
minor, is discussed in relation to graphite intercalated 
compounds (GIG). 

The resulting electronic structure, which agrees well 
with previous density functional theory calculations, [13, -'^^§•[11 



FIG. 2: (Color online) Vibrational properties of picene. From 
top to bottom: Partial C and H contribution to the phonon 
eigenvectors of a picene molecule. Phonon DOS of (undoped) 
solid picene with numbers that mark the energy position of 
selected modes, whose F-point eigenvectors are shown in the 
bottom panel. For clarity, only one picene molecule is shown. 
Modes 4 and 5 give the highest contribution to the ep cou- 
pling, shown in Fig. O 



For all dopings, we will use the phonon spectrum of 
undoped solid picene, whose properties are summarized 
in Fig. [21 Most of the phonon branches of solid picene 
(not shown) have very little dispersion, which reficcts 
the molecular nature of the solid. The dispersive inter- 
molecular modes (acoustic, translations and rotations of 
picene molecules) make up most of the phonon Density 
of States (DOS) for frequencies below ^--^85 cm~^. The 
remaining modes are intramolecular vibration of G and 
H atoms and their general character can be inferred from 
the top two panels of Fig. [2J which show the calculated 
partial eigenvectors of an isolated picene molecule. The 
modes below ^800 cm~^ are due to in- and out-of-plane 
vibrations involving both G and H atoms that cause only 
small changes in the bond length between neighbouring 
atoms. Representative examples arc shown in the bottom 
panel of Fig. [31 They are out-of-plane modes that bend 
picene molecules with frequency of '^55 cm~^ (1) and 
breathing modes with frequencies of ~260 (2) and ~585 
cm^^ (3). At higher frequencies, the vibrational modes 
cause significant change of G-G and G-H bond lengths of 
adjacent atoms. The vibrations between ~800 cm~^ and 
^1300 cm^^ are mostly out-of-plane, while the modes be- 
tween ~^1300 cm^^ and '^1600 cm~^ have G-H bending 
and G-G bending and stretching character. Two typical 
examples of such modes are shown in (4) and (5), with 
frequencies of ^1390 and ~1625 cm~^, respectively. The 
spectral distribution of modes in this region is very simi- 
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lar to that of graphene and (intercalated) graphites. The 
modes very high in energy (~3200 cm~^) have in-plane 
C-H bond stretching character as shown in (6). 

The results of our ep calculations for hole and elec- 
tron doped picene are summarized in Fig. [31 Even using 
the rigid band approximation, obtaining a reliable esti- 
mate of the coupling is far from trivial because the size 
of the system is large. However, the molecular nature of 
solid picene makes it possible to use additional approxi- 
mations, which made the calculations feasible. In normal 
metals, the interaction between phonons and electrons is 
normally expressed in terms of Eliashberg spectral func- 
tion: 
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where 5^ k+q the ep matrix elements that are aver- 
aged on the Fermi surface S{e^). Its first inverse mo- 
ment, A = 2 Jg°° duja'^F{uj)/uj, gives the total ep coupling 
that appears in the exponent of the McMillan formula for 
the critical temperature Tc of ep superconductors. Eval- 
uating expression ( [1]) requires a careful integration in 
reciprocal space over a dense set of phonon (q) and elec- 
tron (k) wavevectors even for elemental metals. How- 
ever, we could obtain well converged results for A and 
a'^F{uj) using affordable q and k meshes by exploiting 
the molecular nature of solid picene. First, we found the 
q dependence of A and a^F^oj) to be weak, reflecting the 
molecular nature of solid picene. Indeed, the ep matrix 
elements g are essentially intramolecular quantities. Ad- 
ditionally, if the Fermi surface is not pathological, the 
integration over k and q in Eq. [1] can be done sepa- 
rately. The dimensionless ep coupling parameter A can 
then be factorized as A = N{Ep)Vep, where the density 
of states N{Ep) is determined by intermolecular interac- 
tions and the ep coupling strength Vgp is essentially an 
intramolecular quantity. [2^, [2^ Similarly, it is possible 
to define a spectral function for the ep coupling strength 
Vep; a'^F{uj) = a'^F{uj)/N{EF)- The two intramolecu- 
lar quantities and Vep converge much faster as a 
function of k-mesh than the intermolecular N{Ep). To 
eliminate the noise coming from pathological regions of 
the DOS, we evaluated Vep and a'^F{ijj) for different fill- 
ings X of the electronic bands as shown in the two insets 
of Fig. [10 

The values of Vep thus determined arc fairly con- 
stant with X within a single subset of bands, as one 
would expect for an intramolecular quantity. We obtain 
Vep = 150 ± 20 meV for holes and Vep = 110 ± 5 meV for 
electrons. This value of coupling is in reasonable agree- 
ment with the empirical trend Vep = 1800/7V^, where iV^r 
is the number of atoms involved in tt states, that has 
been estimated in C tt bonded molecular systems. [2^ 

Before discussing the implications for superconductiv- 
ity, it is insightful to analyze the Vep spectral function 
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FIG. 3: (Color online) Upper Panel: Spectral distribution 
of the ep matrix elements of solid picene, doped with holes 
(green/light gray) or electrons (blue/dark gray). The two 
insets show the doping dependence of the corresponding ep 
coupling strength Vep = / uj~^a'^ F{Lj)dLj for different val- 
ues of doping X. Lower Panel: Electronic DOS per picene 
molecule, in states/(eV spin) as a function of doping x. The 
total ep coupling constant A as a function of doping is given 
by A(x-) = N{x)Vep. 



a^_F(cj), which is shown in the upper panel of Fig. [3] 
The ep spectral function has two main peaks at ^1400 
cm^^ and ~16GG cm^^, which correspond to modes that 
bend C-H bonds and stretch C-C bonds (modes 4 and 
5 in Fig. [2]). Additionally, we also find smaller coupling 
to modes with lower frequencies {uj < 800 cm^^). This 
spectral distribution is in agreement with the model cal- 
culations, (ssj which also predict a sizable contribution at 
^-^260 cm~^ that is present but reduced in our full-linear 
response calculation, [sj] 

The effect of the dopants on ep coupling, which is 
not explicitely included in our calculations, can be es- 
timated from the comparison with graphite intercalation 
compounds. The intercalants in GICs partially reduce 
the ep coupling to C bond stretching modes, but also 
provide some additional coupling to low-lying Einstein 
phonons. 35-33 These two effects balance each other to 
leave A basically unchanged, while win is reduced due to 
the additional weight at low frequencies. In contrast to 
GICs, which have substantial intercalant character for 
states near the Fermi level, both effects should be much 
weaker in solid picene because the metal states are not 
expected to significantly hybridize with C Pz states near 
the Fermi level. [3 ^ 

Accurate values of A = VepN{Ep) as a function of 
doping can be obtained using the DOS shown in the lower 
panel of Fig. [31 which has been obtained using a (8)"^ k- 

For X = 3, we get 
A = 0.78 and win = 
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grid and tetrahedron method 
N{Ef) = 7.06 states spin"! eV 
1021 cm~^ Using ^* = 0.12 in the McMillan formula, we 
obtain a critical temperature Tc ~ 56.5 K. This is more 
than three times the experimental value of Tc, and we 
need to use a larger value of ^* = 0.23 to reproduce the 
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experimental Tc. In the case of hole doping, the T^s arc 
lower due to smaller values of DOS. For x ^ —1.33, with 
the same value of = 0.23, N[Ep) = 4.31 states spin~^ 
eV~^, A = 0.65 and cjin = 890 cm~^, and we obtain a 
maximum Tc of 6 K. 

This relatively high McMillan Tc's that result from 
strong ep coupling to a few high-energy C bond- 
stretching phonons indicates that superconductivity in 
alkali-doped picene and other hydrocarbons is most likely 
phonon-mediated, and related to that of other C and B- 
based superconductors such as MgB2, boron-doped dia- 
mond, and GICs. 

However, the estimates of based on Migdal- 
Eliashberg theory must be taken with care in the case of 
picene. In fact, we have shown that the bandwidth of the 
conduction bands {W ^ 300 meV for electrons and ~ 1 
eV for holes), the frequencies of strongly coupled phonons 
[iOph ^ 200 meV), the ep couphng strength {Vey_ 



110-150 eV) and Coulomb repulsion {U - 1.2 eV|23) 
all have similar magnitudes. In this range of parame- 
ters, the two most important approximations in Migdal- 
Eliashbcrg theory of superconductivity, Migdal's theorem 
and the Morel- Anderson scheme for the screening of the 
Coulomb repulsion, may be invalid. In fuUerenes, whose 
conduction electrons have C Pz character, this regime 
of parameters gives rise to a variety of interesting phe- 
nomena, the most spectacular being the occurrence of 
ep superconductivity near a Mott insulating phase, (soj 
well captured by theoretical studies of local models of 
interacting phonons and electrons in the non-adiabatic 
regime. 0-11 1 

Intercalated picene is also a strong candidate for this 
kind of strongly correlated, non-adiabatic ep supercon- 
ductivity. The experimental situation is still premature 
to verify this claim, but several tests are possible for 
this hypothesis. For example, Migdal-Eliashberg theory 
would predict a smooth behaviour of superconductivity 
as a function of doping and pressure that is governed es- 
sentially by the value of the electronic DOS at the Fermi 
level. On the other hand, strongy correlated local models 
would predict a phase diagram where superconductivity 
exists in close proximity to a possible Mott insulating 
state, which could be tuned by doping or by a change in 
the intermolecular hopping driven by pressure or inter- 
calation of isovalent atoms with different sizes. A very 
exciting prospect is the possibility of studying a wider 
range of intercalated hydrocarbons by changing the num- 
ber and arrangement of benzene rings. This could be an 
interesting avenue for independently tuning the ep cou- 
pling and bandwidth and degeneracy of the conduc- 
tion bands and obtaining an extensive mapping of the 
parameter space for strongly correlated ep superconduc- 
tors. 

In summary, we report the results of our first principles 
calculations for the vibrational spectrum and ep coupling 
of doped solid picene within the rigid band approxima- 



tion. Wc find a large coupling of bond-stretching phonons 
at ~1400 cm~^ and ^1600 cm~^ both to electrons in the 
four lowest conduction bands and holes in the two highest 
valence bands of solid picene {Vep '-^110 and 150 meV, 
respectively). This is sufficiently strong to explain the ex- 
perimental Tc of 18 K reported for K3C22H14 within the 
Migdal Eliashberg theory for ep superconductors. How- 
ever, due to the molecular nature of solid picene, we argue 
that a strongly correlated non-adiabatic ep superconduc- 
tivity similar to that of alkali-doped fuUerides may occur. 
Wc propose possible experimental tests of this hypothe- 
sis and suggest that alkali-intercalated hydrocarbons may 
be ideal playgrounds to study the interplay between elec- 
tronic and vibrational energy scales in a controlled way. 

We would like to thank O. Gunnarsson, R. Arita and 
S. Ciuchi for useful discussion. L.B. acknowledges hos- 
pitality of the KITP, Santa Barbara, where part of this 
work was carried out. 

Recently a Tc of 5 K was reported in K-doped phenan- 
threne (C14H10), confirming the possibility of ob- 

serving superconductivity in other polycyclic aromatic 
hydrocarbons. 

Note added: Wc have become aware of a paper by 
Casula et al. that studies electron-phonon coupling in 
K3C22Hi4.[3H They integrate over full phonon spectrum 
and find A = 0.88 and win = 202 cm^^. With the caveat 
that our low frequency modes are less converged, if we 
also integrate over the full phonon spectrum, we obtain 
A = 1.41 and oj\n = 240.18 cm~^. The discrepancy in the 
value of A may be because our value of DOS (7.06 states 
spin~^ eV~^ per molecule) is more than twice larger than 
theirs (6.2 states spin^^ cV^^ per cell = 3.1 states spin^^ 
cV~^ per molecule). 
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